Dynamics of magnetization on the topological surface 



O 



OS 



Takehito Yokoyama 1 , Jiadong Zang 2 ' 3 , and Naoto Nagaosa 2 ' 4 
1 Department of Physics, Tokyo Institute of Technology, Tokyo 152-8551, Japan 
2 Department of Applied Physics, University of Tokyo, Tokyo 113-8656, Japan 
3 Department of Physics, Fudan University, Shanghai 200433, China 
4 Cross Correlated Materials Research Group (CMRG), ASI, RIKEN, WAKO 351-0198, Japan 

(Dated: March 22, 2010) 

We investigate theoretically the dynamics of magnetization coupled to the surface Dirac fermions 
of a three dimensional topological insulator, by deriving the Landau-Lifshitz-Gilbert (LLG) equation 
in the presence of charge current. Both the inverse spin-Galvanic effect and the Gilbert damping 
coefficient a are related to the two-dimensional diagonal conductivity a xx of the Dirac fermion, 
while the Berry phase of the ferromagnetic moment to the Hall conductivity a xy . The spin transfer 
torque and the so-called /3-terms are shown to be negligibly small. Anomalous behaviors in various 
phenomena including the ferromagnetic resonance are predicted in terms of this LLG equation. 

PACS numbers: 73.43. Nq, 72.25.Dc, 85.75.-d 
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Topological insulator (TI) provides a new state of mat- 
ter topologically distinct from the conventional band in- 
sulator [J. In particular, the edge channels or the surface 
states are described by Dirac fermions and protected by 
the band gap in the bulk states, and backward scatter- 
ing is forbidden by the time-reversal symmetry. From 
the viewpoint of the spintronics, it offers a unique op- 
portunity to pursue novel functions since the relativistic 
spin-orbit interaction plays an essential role there. Actu- 
ally, several proposals have been made such as the quan- 
tized magneto-electric effect [|J,_giant spin rotation Q, 
magneto-transport phenomena [4j, and superconducting 
proximity effect including Major ana fermions 0- 



Also, a recent study focuses on the inverse spin- 
Galvanic effect in a Tl/ferromagnet interface, predicting 
the current-induced magnetization reversal due to the 
Hall current on the TI [8( . In Ref . [8j , the Fermi energy is 
assumed to be in the gap of the Dirac dispersion opened 
by the exchange coupling. In this case, the quantized 
Hall liquid is realized, and there occurs no dissipation 
coming from the surface Dirac fermions. 

However, in realistic systems, it is rather difficult to 
tune the Fermi energy in the gap since the proximity- 
induced exchange field is expected to be around 5-50meV. 
Therefore, it is important to consider the generic case 
where the Fermi energy is at the finite density of states 
of Dirac fermions, where the diagonal conductivity is 
much larger than the transverse one, and the damping of 
the magnetization becomes appreciable. Related systems 
are semiconductors and metals with Rashba spin-orbit 
interaction, where the spin-Galvanic effect and current 
induced magnetization reversal have been predicted Q 
and experimentally observed [Tot Compared with 

these systems where the Rashba coupling constant is a 
key parameter, the spin and momentum in TI is tightly 
related to each other corresponding to the strong cou- 
pling limit of spin-orbit interaction, and hence the gigan- 
tic spin-Galvanic effect is expected. 
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FIG. 1: (Color online) (a) Illustration of the Dirac dispersion 
on top of TI. The Fermi level £f is far above the surface 
gap opened by magnetization in the ferromagnetic layer, (b) 
Sketch of FMR experiment in the soft magnetic layer. The 
substrate in the figure is TI, which is capped by a layer of 
soft ferromagnet. The magnetization precesses around the 
external magnetic field H e g. 



In this letter, we study the dynamics of the magnetiza- 
tion coupled to the surface Dirac fermion of TI. Landau- 
Lifshitz-Gilbert (LLG) equation in the presence of charge 
current is derived microscopically, and (i) inverse spin- 
Galvanic effect, (ii) Gilbert damping coefficient a, (iii) 
the so-called /3-terms, and (iv) the correction to the Berry 
phase, are derived in a unified fashion. It is found that 
these are expressed by relatively small number of param- 
eters, i.e., the velocity vf, Fermi wave number fcj?, ex- 
change coupling M , and the transport lifetime r of the 
Dirac fermions. It is also clarified that the terms re- 
lated to the spatial gradient are negligibly small when 
the surface state is a good metal. With this LLG equa- 
tion, we propose a ferromagnetic resonance (FMR) ex- 
periment, where modifications of the resonance frequency 
and Gilbert damping are predicted. Combined with the 
transport measurement of the Hall conductivity, FMR 
provide several tests of our theory. 
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Derivation of LLG equation. — By attaching a ferro- 
magnet on the TI as shown in Fig. [1] we can consider a 
topological surface state where conducting electrons in- 
teract with localized spins, S, through the exchange field 



M J dr n(r) • cr (r) . 



(1) 



Here, we set S = Sn with a unit vector n pointing in the 
direction of spin, cr(r) = c' (r)crc(r) represents (twice) 
the electron spin density, with c 1 = (c|, ct) being electron 
creation operators, cr the Pauli spin-matrix vector, and 
M being the exchange coupling energy. The total Hamil- 
tonian of the system is given by H to t = H$ + H e \ + H ex , 
where H$ and H e \ are those for localized spins and con- 
ducting electrons, respectively. 

The dynamics of magnetization can be described by 
the LLG equation 



n = 7oHeff x n + a^h x n + t' cl 



(2) 



where 7off e ff and cto are an effective field and a Gilbert 
damping constant, respectively, both coming from H$. 
Effects of conducting electrons are contained in the spin 
torque 



t e l(r) = s t' el (r) = Mn(r) x (cr(r)) r 



(3) 



which arises from H ex . Here, sq = S/a 2 is the local- 
ized spin per area a 2 . In the following, we thus calculate 
spin polarization of conducting electrons perpendicular 
to n, (&±(r)) ne , in such nonequilibrium states with cur- 
rent flow and spatially varying magnetization to derive 
the /3-term, or with time-dependent magnetization for 
Gilbert damping. Here and hereafter, (• • • } nc represents 
statistical average in such nonequilibrium states. 

Following Refs. flil - [l4j we consider a small transverse 
fluctuation, u = (u x ,u y ,0), \u\ <C 1, around a uniformly 
magnetized state, n = z, such that n = z + u. In the 
'unperturbed' state, n = z, the electrons are described 
by the Hamiltonian 



Ho = {kyO- 1 



Ma z -e F + V imp (4) 



where Vi mp is the impurity potential given by Vi mp = 
u^ i 5{r — Ri) in the first-quantization form. We take 
a quenched average for the impurity positions Ri. The 
electron damping rate is then given by 7 = l/(2r) = 
imiU 2 VF in the first Born approximation. Here, rij is the 
concentration of impurities, and v F = £p / (2irv 2 F ) is the 
density of states at £f- We assume that 7 vp^F = 
\Js 2 F — M 2 , M, and calculate spin transfer torque in the 
lowest non-trivial order. 

In the presence of u(r,t) = u(q,uj) e Ki- r - u t) } the con- 
ducting electrons feel a perturbation (note that H e \ + 
Hex = Ho + Hi) 



Hi 



■u(q,w) 



(5) 



and acquires a transverse component 

(a'Z{q,u)) ue = Mx¥{q,w + to)u p {q,u) (6) 

in the first order in u in the momentum and frequency 
representation. Here, x± ls the transverse spin suscep- 
tibility in a uniformly magnetized state with a,0 = x, y, 
and summing over (3 is implied. 

Now, we study the w-linear terms in the uniform (q = 
0) part of the transverse spin susceptibility, xl^iQ = 
0, u) + iO). We make the following transformation of the 
operator: 



c= Uc = 



v /2e(e + M) 



v F (k y + ik x ) 
e + M 



(7) 



with e = y/(v F k) 2 + M 2 . Note U^U = l,U^a x U = 
vpky/e, and U'o~ v U = —v F k x le. This transformation 
maps two component operator c into one component op- 
erator on the upper Dirac cone c. With this new op- 
erator, we calculate the transverse spin susceptibility in 
Matsubara form 



X± (0, jwa) = 



/ dre l ^ r (ryfo.Vfo.o)) 

Jo 



= -T 



U^a a UG(k, ie n + iuj x )U^a p UG(k, ie n ) (8) 
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with G(k, ie n ) = (ie n - e + e F + i7Sgn(e„)) _1 . By sym- 
metry consideration of the integrand in fe-integral, we 
find xl /3 (0 , ico\) cx 6 a fj. After some calculations, we ob- 
tain the torque stemming from the time evolution: 



, iuj I / v F kf 



2n 2v F \ e F 
_ 1 / Mv F k F 
= 2 \~Tf~ 



epTri x u 



vpTW x n. 



(9) 
(10) 



This result fits the conventional Gilbert damping with 

2 _2 



1 / Mvpkr 
a = - 



2 V £F 



VpT 



ns' 



(11) 



We next examine the case of finite current by applying 
a d.c. electric field E, and calculate a linear response of 
o°t to E, i.e., < crf(q) > ne = Kf{q)E l . First, it is clear 
that Kf(q = 0) = -e ia a xx /(ev F ) where e ia and a xx 
are the anti-symmetric tensor and diagonal conductivity, 
respectively, because electron's spin is "attached" to its 
momentum. This represents the inverse spin-Galvanic ef- 
fect, i.e., charge current induces magnetic moment. Since 
we assume that Fermi level is far away from the surface 
gap, o~ xx 3> o~ X y where a xy is the Hall conductivity. The 
dominant term in x is thus Xxy oc o~ xx . This is quite 
different from the case studied in Ref. [8j, where Fermi 
level lies inside the surface gap and therefore a xx is van- 
ishing. Hence, the only contribution to the inverse spin- 
Galvanic effect is Xxx oc cr xy , which is much smaller than 
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the effect proposed in this letter. Compared with the in- 
verse spin-Galvanic effect in Rashba system (§-l_l|, this 
effect is much stronger since the small Rashba coupling 
constant, i.e., the small factor a^kp/Ep in Eq. (16) of 
Ref. 0, does not appear in the present case. Taking into 
account the realistic numbers with a = lO~ 11 eVm and 



vp = 3 x 10 5 m/s, one finds that the inverse spin-Galvanic 
effect in the present system is ~ 50 times larger than that 
in Rashba systems. 

The next leading order terms of the expansion in u@ 
and Qj ca n be obtained by considering the four-point ver- 
tices [12j as 



J 



7r hi 



4 %-KEp 



£ik£ 3 l [Sasaki + SakSpi + S a l5/3k] qjU^Ei 



-eM 



hi 
324 



[q • Eu" - q • (u x z)(E x z) a + u • (E x z)(q x z) a ] 



Therefore, the spin torque steming from the spatial gradient has the form: 

if, = -p— [n x (j ■ V)n - (j - (j • n)l)V • (n X z) + ( V - (n • V)z)n • (j x z)] 



(12) 
(13) 

(14) 



where j = er^E with charge current j and conductivity 



e 2 ( vpk F \ ' 
4ir \ e F J 



SpT. and 



5?r 



M 



4epT \vpk F 



(15) 



From Eq. ([M)) . one can find the followings: (i) The spin 
transfer torque of the form (j • V)n is missing since 
we consider the upper Dirac cone only, (ii) The (3- 
tcrm has a form essentially different from that in the 
conventioal one. 12 , [Til Il6| In contrast to the conven- 
tional ferromagnet. [lathis constant comes from the non- 
magnetic impurity. Considering vpk F = ep, we get 
a/P ~ (sft) 2 from Eqs. (JTTJ and ([T5]l. Therefore, 
the /3-terms are negligible for a good surface metal, i.e., 

SpT > 1. 

Up to now, we consider only one branch of the band 
where the Fermi energy is sitting. When we consider the 
2-band structure, i.e., the 2x2 matrix Hamiltonian TL = 
v F [(k y + ^^)cF x - (k x - MJhL^y] ) we have the correction 



to the Berry phase term. In analogy with the minimal 
coupling of electromagnetic field, A = — j^(— n y ,n x ) 
plays the same role as the U(l) gauge. By integrating 
the fermions out, one can get a Chern- Simons term in 
terms of the magnetization Lcs = &xy^ vp A^dvAp where 
fx, v, p = t, x, y. When the gradient of magnetization van- 
ishes, it can be rewritten as 



Lcs = <rx V {—) 2 {n x h y - n y h x ). (16) 

tVp 



This additional term can be interpreted as an additional 
Berry phase for the magnetization. In fact, as n z remains 
constant in the present case, we have [n x , n y ] = in z . 
Therefore, n x and n y become conjugate variables up 
to a factor, which naturally leads to a Berry phase: 
n x h y — n y h x . This term is exactly equivalent to the 
Chern-Simons term. 

Including all the terms derived above, we finally arrive 
at a modified LLG equation: 



n - 2a xy {— ) 2 n/(s N) = j H cff x n - 
evp 



M 



evFSoN 



(-j + (n ■ j)z) + (ao + a/N)n xn + t p J(s Q N) (17) 



r 



where N is the number of ferromagnetic layers. Note that 
a-, ft- and Berry phase terms originate from the interplay 
between Dirac fermions and local magnetization which 
persists over a few layers of the ferromagnet. Therefore, 
the overall coefficients are divided by the number of fer- 
romagnetic layers N. 



Ferromagnetic resonance. — Observing the small value 
of /3, the spatial gradient of magnetization can be ne- 
glected for the time being. Only one uniform domain in 
the absence of current is taken into account for simplic- 
ity. Without loss of generality, assume that an external 
magnetic field is applied along z direction, and consider 
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the ferromagnet precession around that field. n z = 
is kept in the first order approximation, namely n z is a 
constant in the time evolution. By inserting the ansatz 
n x(y){t) = n x(y) e ~ lut m to the modified LLG equation, 
one obtains 



e + n 



(18) 



where r\ — (ao + a/N), ujq — 7o-ff e // and £ = 1 — 
2a xy{-^;) 2 1 '(sqN). Expanding up to the first order in 
a xy and rj, one gets Sftw = lu q + 2a xy (-^) 2 uj Q /(s N) 
and So; = rjuiQ. Therefore, the precession frequency ac- 
quires a shift proportional to a xy in the presence of in- 
terplay between Dirac fermions and the ferromagnetic 
layer. The relative shift of is 2a xy (-^j-) 2 Ld /(soN) = 

dwi7(^r) 2 ~ Ml7 ) 3 0- By tuning the Fermi level, 
this shift can be accessible experimentally. 

Meanwhile, the Gilbert damping constant a can be 
measured directly without referring to the theoretical 
expression in Eq. (jlljl - One can investigate the fer- 
romagnetic layer thickness dependence of FMR line- 
width. While increasing the thickness TV of ferromagnet, 
the Gilbert damping constant stemming from the Dirac 
fermions decreases inversely proportional to the thick- 
ness. Taking into account the realistic estimation with 
eft ~ 100 and M/ef ~ 0.3, one has a/so ~ 1, while 
ao ~ 0.001 usually. Therefore, even for a hundred of lay- 
ers of ferromagnet, the contribution from the proximity 
effect is still significant compared to the one coming from 
the ferromagnet itself. Observing that the imaginary part 
of resonance frequency in Eq. (1181) is proportional to 77, 
one may plot the relation between the FMR peak broad- 
ening, namely SScu, and 1/N. The broadening is a linear 
function of 1/N, and approaches the value of the ferro- 
magnet at large thickness limit. We can find the value of 
a from the slope of the plot. 

On the other hand, the real part of FMR frequency 
provides rich physics as well. Since in the presence of 
additional Berry phase, the frequency shift is propor- 
tional to the Hall conductivity on the surface of TI, it 
leads to a new method to measure the Hall conductiv- 
ity without four-terminal probe. In an ideal case when 
the Fermi level lies inside the surface gap, this quantity 
is quantized as a xy = However, in realistic case, 

Fermi level is away from the surface gap, and therefore 



the Hall conductivity is reduced to a xy 



Mr. 



ll 



As 



W ~ 2h e F 

a result, the shift of resonance frequency is proportional 
to n 2 oc cos 2 9, and the FMR isotropy is broken. Here, 
9 is the angle between effective magnetic field and the 
normal to the surface of TI. One can perform an angle 
resolved FMR measurement. The signal proportional to 
cos 2 9 comes from additional Berry phase. 

Since parameters a and j3 depend on M and r, it 
is quite important to measure these quantities directly. 



Molecular-beam epitaxy method can be applied to grow 
TI coated by a thin layer of soft ferromagnet. As is re- 
quired in the above calculation, Fermi level of TI should 
lie inside the bulk band gap. Also, the soft ferromag- 
net should be an insulator or a metal with proper work 
function. One may employ angular resolved photoemis- 
sion spectroscopy (ARPES) or scanning tunneling micro- 
scope techniques to measure the surface gap A opened 
by the ferromagnet, which is given by A = Mn z . As the 
easy axis n z can be found experimentally, M can be fixed 
as well. On the other hand, the lifetime r is indirectly 
determined by measuring the diagonal conductivity a xx 

via a xx — ( " J £ft. Finally, Fermi surface can 
be determined by ARPES, and all parameters in LLG 
equation Eq. fTTf can be obtained. 

In summary, we have investigated theoretically the 
dynamics of magnetization on the surface of a three 
dimensional topological insulator. We have derived 
the Landau-Lifshitz-Gilbert equation in the presence of 
charge current, and analyzed the inverse spin-Galvanic 
effect and ferromegnetic resonance predicting anomalous 
features of these phenomena. 
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